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After revisiting the dynamical system based approach of spherically symmetric vacuum 
braneworlds, pointing out and studying the existence of a transcritical bifurcation as the 
dark pressure parameter changes its sign, we analyze some consequences of not discard the 
brane cosmological constant. For instance, it is pointed that the existence of an isothermal 
state equation between the dark fluid parameters cannot be obtained via the requirement of 
a quasi-homologous symmetry of the vacuum. 

I. INTRODUCTION 

Since 1999, braneworld models has attracted much attention of the scientific community. In fact, 
the view of our universe as a four-dimensional brane embedded into a five dimensional warped bulk 
[T] leads to new insights concerning well posed problems from particle physics to cosmology [2]. 
Several aspects of the braneworld cosmology were scrutinized, such as the gravitational collapse 
[3], the cosmological dynamics [3], inflation [5j and cosmological perturbations [6j (for a broader 
review, see [?]). A quite interesting approach, opening several new perspectives was developed in 
Refs. [51 [9]. In fact, these works lead to an effective braneworld gravitational equation encoding 
local as well as non-local (purely geometrical) corrections. The non-local corrections are given by 
the explicit appearance of a specific bulk Weyl tensor projection. 

In order to study the non-local contribution of the resulting gravitational equation, it is usual 
to parameterise the Weyl tensor by a type of (dark) cosmological fluid, respecting all the necessary 

* Also at The Niels Bohr Institute - Copen hagen University, Blegdamsvej 17, DK-2100 Copenhagen, Denmark. 
'Electronic address: mabdalla@ift.unesp.br 

^Electronic address: pablofisico@iftTunesp.br' 

^Electronic address: hoff@feg.unesp.br; hoff@ift.unesp.br 



2 



constraints of the projection procedure. An analysis of spherically symmetric branes in the vacuum 
state, specially devoted to some interesting brane symmetries investigation, was performed in [10] . 
In particular, in this last paper, the notion of quasi-homologous branes, i. e., branes invariant 
with respect to the group of quasi-homologous transformations, is presented. In Refs. |1H - I13| . a 
qualitative analysis was carried out, based upon the dynamical system defined by the projected 
equations. 

In this paper we firstly aim to investigate some qualitative aspects of quasi-homologous spher- 
ically symmetric vacuum braneworlds without an effective cosmological constant. To some extent 
it was previously developed [lOj . nevertheless we shall revisit this program, calling attention to 
the existence of a transcritical bifurcation as the dark pressure parameter is swept. This analysis 
is performed in Section III, being Section II devoted to some review and first steps towards the 
structure equations. Then we move forward, rewriting the vacuum structure equations, this time 
with a brane cosmological constant. It is shown in Section IV that, in this last case, the effective 
cosmological constant, although quite small, may lead to important consequences in the braneworld 
picture. In particular, it is not possible to ensure an isothermal equation of state for the dark fluid 
by requiring a quasi-homologous invariant vacuum. Two aspects, in this research line, shall be 
mentioned. Firstly, an argument in favor of keeping the cosmological constant term in the analysis. 
The usual approach is to discard the contributions coming from A. In fact, keeping in mind the 
minuteness of the brane (4D) cosmological constant, its effects can be safely neglected in finite 
sized gravitational systems as stars, galaxies, etc. However, since we are dealing with the vacuum 
of the brane itself, to give up of the brane cosmological constant terms seems to be an oversimpli- 
fication. We note, in advance, the appearance of several terms involving a product of A with the 



radial coordinate in the formulae (10)-(13) below, making expected some effect in regions where 
r is very large. Secondly, the fact that an isothermal equation of state for a cosmological fluid is 
consequence of a quasi-homologous symmetry (for usual systems) was reported long ago [H] and 
translated to the braneworld language more recently p^. Hence, this type of symmetry breaking 
due to the presence of the cosmological constant is an important improvement in the braneworld 
scenario study. 



II. THE STRUCTURE EQUATIONS 



It is well known that the projection of the 5D gravitational field equations on the brane leads to 
important corrections classified in two classes: the local corrections, encoded in the quadratic brane 
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stress tensor terms, and the non-local (purely geometric) ones, encrypted in a specific projection 
of the bulk Weyl tensor. More precisely, being the five dimensional Einstein equations given by 

Gij = kITij, (1) 

where Gjj is the Einstein tensor, Kg = SttG^, and Tjj = A^gjj + 6{y)[—Xhgij + t/j], it is possible 
to show that the Gauss-Codazzi procedure leads to the following gravitational equations on an y 
constant hypersurface [8] 

Gfj,u = —Ag^iu + i^\'TiJ,v + i^t^iJ.i' ~ ^iiv- (2) 

The notation is now evident: capital Latin indices go from 0..4, while Greek indices vary in the 
range 0..3. In the above expression, S^y is the aforementioned tensor encompassing quadratic brane 
stress tensor terms, as our analysis is in the brane vacuum this tensor vanishes. Before exploring 
the Weyl tensor E^y however, it is important to remark that after the projection procedure the 

2 

constants A and K5 are no longer fundamental. Instead, they are given by A = -^(As + k|A^/6) 
and — Ave Ab/6, where A;, is the brane tension. In the vacuum, the Eq. ([2]) reduces to 

R^lv = —Efj,y + Agf^y. (3) 

As the standard approach to the Weyl tensor, its is parameterized in the following way |15j 
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(4) 
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This specific form for Ef^y is motived by the symmetries that it must obey. In Eq. Q li'^ is a given 
four velocity field and /i^;^ = gf^y + Ufj,Uu. It is possible to show [lU] that for a static spherically 
symmetric vacuum on the brane, Q/^ = and P^y = P{r){r^ry — /i^i^/3), where is the unit radial 
vector and r the radial distance. Besides, in an inertial frame on the brane one has = (1,0) 
and h^i, = diag{0, 1, 1, 1). In the current jargon the decomposition Q is the so-called Weyl fluid, 
U = U{r) is the dark radiation, and P{r) is said as the dark pressure. This function will play an 
important role in the qualitative analysis to be presented in the next section. 

After these introductory review we shall reobtain the structure equations for this case without to 
discard the brane cosmological constant. Using the usual static spherically symmetric line element 
on the brane 

ds^ = -e^^^^dt^ + e^^'Ur'^ + r\dd^ + sin'^ed^'^), (5) 

we have 
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e'^l \ + 4 = 3aU + A, 



(6) 



u' 1 \ 1 

7 + ? " 



a(?7 + 2P) - A, 



(7) 



q(C/- P) - A, 



(8) 



6P 



2C/ + P r(2[/ + P)^ 



(9) 



where a = ^^4^^"^ and a prime denotes derivation with respect to r. From Eq. (|6|) one can see that 



e-A = 1 _ C _ Q _ AH 
r r 3 



being Q = 3a / Ur^dr. Substituting Eq. ^91) in Q one has 



ap 1 

\J' = -2P' - — + - (1 - e^)(2C/ + P) - (2C/ + P)re^[a(f/ + 2P) - A] 



(10) 



(11) 



and, with the aid of Eq. (|10p, one arrives at the following expression 
dU 



p {2U + P){C + Q + ar-'{U + 2P)-lAr^ 

2— -6 

dr r 



dr 



r2 1 



c 



(12) 



This last equation, together with ^ = 3ar^f7, will be extensively studied in the following Sections. 



III. QUASI-HOMOLOGOUS BRANES: DYNAMICAL SYSTEM ANALYSIS 



In order to set up a dynamical system analysis of the vacuum equations let us define 

C Q Ar2 

q = - + - + ^, 

r r 6 



(13) 
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together with the change of coordinate 6 = ln{r). It can be readily verified that 

dq 



while a bit of algebra leads to 



- = 2M-2p-2--12e A — . (15) 



We note that, due to the presence of the cosmological constant term, the equations (14) and (15) 
are slightly different from the previously ones obtained in the literature. 

After all the comprehensive discussion presented before |10| I13j . we would like to revisit the 



qualitative analysis of Eqs. (14) and (15), discussing another way to reach a quasi-homologous 
system, as well as interpreting the behavior of the fixed points. In order to accomplish that, 
we shall disregard the cosmological constant term along the present Section. As mentioned in the 
Introduction, when dealing with the vacuum of the brane itself, it seems to be an oversimplification. 
Even so, one may acquire physical insight by performing such a simplified analysis. We postpone 
to the next Section the importance of the cosmological constant for such a system. 



Without the brane cosmological constant term, Eqs. (14) and (15) provide an useful starting 
point for the dynamical system analysis. Note, however, that the dp/dO term is an obstruction to the 
dynamical system program (it is well known that dynamical systems in three or more dimensions 
may behave exotic. For instance, it was demonstrated in [16] that the generic behavior of three 
dimensional trajectories approachs some strange attractor, although all the possible behaviors were 
not mapped yet). Generally, the approach used to overcome this problem is to implement ah initio 
an equation of state between U and P (^u and p) . This procedure is closely related to the existence 
of a quasi-homologous symmetry on the brane. In fact, in Ref. [10] it was stated that a theorem 
asserting that the vacuum brane equations (in the case we are dealing with, i. e., spherically 
symmetric and static) are invariant with respect to the group of quasi-homologous transformations 
if and only if P = 7{7, where 7 is a constant. We shall comment more on that theorem in the 
next Section. Our procedure now is, instead of implementing the quasi-homologous symmetry as 
an input, to require a constant p, since in this way we immediately get an autonomous dynamical 



system out of (14) and (15). Obviously, it is also an assumption, and at first sight it is not related 
to the usual one {P = 7C/) which engenders the quasi-homologous symmetry. It is possible to see, 
however, that the underlying dynamical system leads to the case P = 7C/. 

In order to accomplish that, we rewrite the structure equations with all the simplifications taken 
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into account. Eq. (14) remains unchanged, where now q = ^^'^ and /i = Sar'^U, and 

i2f,+p)[{fi + 2p)/3 + q] 



de 



2fi-2p 



1 



(16) 



being p = 3ar^P. Notice that the assumption that p is constant means P ~ Let us write 



p = 3a/3 with /3 constant {P = fi/r"^)- Eq. (14) trivially gives a strong constraint which must be 
respected by every single fixed point in the (g, /u)— plane, namely /i* = (?* where * labels a fixed 



point. The fixed points are completely obtained from Eq. (16). Hence, it reads 

3 
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We note from (17) that 



a/3(9a/3 + 10) < 1/3, 



(17) 



(18) 



in order to have a real fixed point. This constraint together with the fact that /z* 7^ 1 (see Eq. 



(16)), will play an important role in the subsequent analysis. 
Notice that from /i* = we have 



Nevertheless, by the definition of Q, we have 



3ar2[/ 



dQ 

dr ' 



in such a way that 



M = C + Q. 
dr 



(19) 



(20) 



(21) 



Taking the derivative on both sides of Eq. (21) with respect to r, we get 



dr"^ 



0, 



resulting in 



dQ 
dr 



Sar^U = (5, 



where 5 is a constant. As /x* = (7* for P = ^ {p constant) it is readly verified that 



(22) 



(23) 



which means the existence of the quasi- holonomic symmetry jlOj . 
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Now, following the standard classification of critical points [T7], it is possible to find that the 
eigenvalues of the Jacobian matrix is given by 



1 



19^ + ^ - ^ ± ^193/i2 + 570/u*a/3 - 6/i, + Ula'^p^ + 18a/3 + 9 



(24) 



Hence it is necessary to look at the real part of A for all the cases, namely {fif,X-^-), {fj,f,X-), 
and (//~,A+), (;U~,A_). The classification is quite characteristic [17]: in a given fixed point, if 
Re{X±) > 0, then the fixed point is a repellor; if Re{X±) < 0, then the fixed point in question is 
an attractor. Finally, if i?e(A-|_) and i?e(A_) have opposite signs, then the fixed point is called a 
saddle point. As mentioned, we shall make explicit the existence of a transcritical bifurcation in 
the dynamical system in question. In order to unveil this characteristic, we shall investigate the 
behavior (classification) of each critical point as the dark pressure parameter changes its sign. As 
we shall see, there is a complete change in the behavior of one of the critical points, characterizing 
the aforementioned bifurcation. 



A. The p = case 



In order to get a physical insight on the possibilities raised in the scope of the dynamical systems 



behind the structure equations, let us investigate the simplest case with p = (/? = 0) in Eqs. (17) 



and (24). From (17) it is simple to see that the critical points are fif = 3/7 and /i^ =0. For the 



case /i+ = 3/7, the eigenvalues (24) are given by 

1 



and 



A+ = --(-2 + V228) < 
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A_ = --(-2 - V228) > 0, 



therefore this point is a (less important) saddle point. The other critical point (^^ = 0) has 



A- 



2A. 



2, 



(25) 



(26) 



(27) 



hence it is a repellor. Note that from the constraint /i* = = p = we need necessarily to have 
C = 0. This fact is indeed technically sound, since in the absence of the dark fiuid we should expect 
a complete vacuum on the brane, and C = (when comparing with the Schwarzschild case) shall 
be recognized as 2MG = 0. In this specific critical point, the gravitational equations give 



svn 



(28) 
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with a constant absorbed in the t coordinate. This "spherical symmetry", from the analysis of the 
dynamical system, is obviously unstable. 



B. The p> case 



The fisrt thing we shall note in this case is the constraint (18). Since now /3 > 0, the allowed 



range of values oi a/3 is within the interval ( 



-30+v/1008 
54 



stressing for the minuteness of possible 



values. Therefore, it is necessary to substitute the Eq. (17) into (24) analyzing all the possible 



cases. This approach is, however, a hard task. Hence it is useful to plot all the possible cases and, 
then, scrutinize the sign possible variations of the eigenvalues. 

From the Figs. (1) and (2), it is possible to conclude that for the fixed point fif the eigenvalue 
A_ is always positive within the relevant range. Instead, A„ can be positive or negative. Therefore, 
fif is either saddle or repellor. 

The situation is also simple for the fixed point. From Figs. (3) and (4) it can be readly verified 
that this point is always a repellor. Nevertheless, the situation can be modified if we consider the 
negative reduced dark pressure case. 




FIG. 1: Eigenvalue A+ for the fixed point /i^. 
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FIG. 4: Eigenvalue A_ for the fixed point fi^ . 



C. The p <0 case 



The last possibility in the possible range of the reduced dark pressure is the negative case. In 



this case, the constraint (18) reads 



a|^|(10-9a|;g|) > (29) 



30+V1008 



Notice the peculiar 



thus we can see that the relevant range for this case is a\f3\ £ |^0, 
behavior presented in Figs. (5) and (6). Inspite of the absence of similarity between the curves, 
there is an interesting pattern around the value 0.7 of the horizontal axis. The change of sign occurs 
at the same interval, in such a way that the critical point /x^ is always a saddle point. Another 
remarkable behavior is happening, this time with the point fi^ . The Figs. (7) and (8), show that the 
eigenvalues at this last critical point are both positive until p/3 = 1, and therefore it is classified as 
a reppelor. Notwithstanding, from this point on, both eigenvalues become negative, and the critical 
point /i~ starts to act as an attractor. Comparing with the p > case, in which this fixed point 
is always a reppellor, now we have a drastic change of classification. This typical behavior is well 
known in dynamical system theory, and characterizes the so-called transcritical bifurcation. 
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FIG. 7: Eigenvalue A+ for the fixed point . In this plot, the horizontal axis has been shortened for the 
sake of exposition. 
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FIG. 8: Eigenvalue A_ for the fixed point /i,^ . In this plot, the horizontal axis has been shortened for the 
sake of exposition. 
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IV. GENERAL SYSTEM WITH BRANE COSMOLOGICAL CONSTANT: BREAKING 

USUAL QUASI-HOMOLOGOUS SYMMETRY 



It was first noted by C. B. Collins |14j that the arguments, somewhat ad hoc, leading to the 
equation of state of static (Newtonian and relativistic) stars may be obtained from the requirement 
that a (quasi) homologous family of solutions should exist for this problem. This concept was 
successfully extended to the (spherically symmetric and static) braneworld case in [lOJ. In this 
last case it was shown that the infinitesimal operator generating the group of quasi-homologous 
transformations is given by 

« = '^-^^ilj + <«+^)<|- (^°' 

Moreover, it is shown that the invariance with respect to this quasi-homologous transformations 
ensures the relation P = 'jU, with 7 constant, and vice- versa [lOj . 

We shall demonstrate here that the usual quasi-homologous symmetry is breaking due to the 
presence of the cosmological constant. Hence, in particular, the simple equation of state is not more 
sufficient to ensure the quasi-homologous symmetry. Before doing so, however, let us make a brief 
account on the aforementioned symmetry. 

Generally speaking, it is said that a given system of (ordinary) differential equations 

Ay™ 

^ = r(t,u), (31) 

being m = 1, ..,n and u = (n^, ..,«") is invariant under quasi-homologous transformations if, and 
only if 

d7/"^(u™) dTr{t) 



du^ dt 



H(/n|/-|), (32) 



with no sum in the m index. In the above expression, H = 7r(t, u)-^+rj"^{t, u) is the infinitesimal 



generator whose action leaves the set (31) invariant. The reason for the specific form of Eq. (32) 
rests upon the application of Lie group technics in the investigation of differential equations systems 
and we refer the reader for more details, for instance, to the Ref. [18]. It is worth, however, giving a 
precise definition for quasi-homologous symmetry: in the above context it means symmetry under 
transformations as t — t- t{t) and u*" — t- u"^{u"^). After all, the constraints imposed by the left 



hand side of Eq. ( 32 ) on the allowable quasi- homologous transformations are such that the unique 
possibility is given by a simple rescaling [14J. This fact will be important in the interpretation of 
quasi-homologous symmetry breaking. 
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Returning to our problem, let us assume that the dark pressure and the dark energy are related 
by the simple state equation P(U) = 'jU, being 7 constant. Remember that, as mentioned, it 
is a sufficient condition to ensure quasi-homologous symmetry in the case without cosmological 
constant. The relevant vacuum gravitational equations are given by 



dQ 
dr 



(33) 



and 



dU 
dr 




[1 + {C + Q + ar^Uil + 27) 



(34) 



Let us investigate the possibility of the usual quasi-homologous transformations generated by the, 
also usual, infinitesimal generator 



n = e(r)^+v'(U)^,+v'(Q)^, (35) 
deliberately taken as being the same of Refs. [IOl[llj (it is in this sense we call it 'usual'). Applying 



the generator (35) in (33) we have, from (^32|, the result 



drj'^iQ) _ dew ^ 2^(r) , rj^U) 
dQ dr 



+ 



U 



(36) 



whose solution is given by [TO 



e(r) 



+ br, 



(37) 



ryi([/) = {c-3b)U 



(38) 



and 



r]'^iQ) = cQ + d, 



where a and d are integration constants, whereas c and 36 are separation constants. 



(39) 



The equation to be satisfied in the case of Eq. (34) is more involved. The left hand side, given 



bv 



"^^jp-, is trivially obtained from Eqs. (37-39). It reads simply 



2a 



+ c-Ab. 



(40) 
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The complete Eq. (32) in this case reads 

(1 + 27-i)(l + 2-f)3ar^U + 6 - 4Ar2(2 + 



2a 



+ c - 46 



+ br 



1+1 



(C+Q1\ 4Ar2 



Fa 



+ {c-3b)U 



1 ar3(l + 27~i)(l + 27) 



+ (cQ + d) 



27-1 - 5 



+ 



(41) 



where Fa = (1 + 27-i)[C + Q + ar3[7(l + 27)] + 6r - 6(C + Q) - ^(2 + 7-I). After a simple, 



but lengthy, algebra it is possible to rewrite Eq. ( 14 ) as 

C Q 



1 



C Q 



3a 



1 



2a , \ . 9 f d cQ \ a 
-. + b]Ar^- - + — +^ + 6 



3a 



+ cr^ (1 + 27-i)(l + 2j)ar'^U + 6(^^ + + (cQ + d){2-f-^ - 5) - 4Ar2(2 + 7-^) (^^ + br^ 



(1 + 27-i)(l + 27)ar3C/ + 6r + (C7 + Q)(27-i - 5) - ^(2 + 7-I) 

It can be easily verified that, as expected, iia = b = c = d = the above constraint is satisfied. 
In fact, it corresponds to the identity transformation (H = 0). Another important remark is that 
if we take A = 0, then the only possibility to satisfy the constraint is a = 0, d = C leading to 
c = b = 1, which is exactly the solution found in Ref. [inj, linking, in this way, the state equation 



with the generator (30). Note, however, that the presence of the brane cosmological constant 



term prevents the existence of (integration and separation) constants that satisfy the equality, 
breaking the quasi-homologous symmetry. We shall interpret this fact as follows: as mentioned, the 
existence of a quasi-homologous symmetry means invariance under simple rescaling of the physical 
parameters. Nevertheless, the cosmological constant has a small fixed value. In the braneworld 
paradigm it depends on the brane tension and on the five-dimensional cosmological constant, both 
fixed parameters. Besides, this dependence is such that it enables a quite specific value for the 
brane cosmological constant. Therefore, from this perspective, the system of equations taking into 
account the cosmological constant term shall not be invariant under rescaling transformations. 



V. CONCLUDING REMARKS 



We would like to emphasize in this concluding remarks the main accomplishes of this work. By 
investigating the subjacent dynamical system associated to the brane vacuum equations for the 
spherically symmetric and static case, we called the attention to the existence of a transcritical 



16 



bifurcation by the variation of the dark pressure parameter. It is important to stress that, in view 
of Section III, the dark pressure parameter must be negative (the existence of the attractor point 
occurs only in this case). Besides, it was shown the possible values oi p by inspecting the necessary 



constraint ( 18 ) to have a real fixed point in the configuration plane. It is important to call attention 
to the minuteness of p. Note also that in all the cases we have 7^ 1 for the allowed p values, which 
is desirable for inner consistency of the dynamical system. Furthermore, we shown that the usual 
equation of state P = (or the existence of quasi-homologous symmetry) is also a consequence 
of the requirement p constant. 

Going further, it was demonstrated that the presence of the brane effective cosmological constant 
term (A) breaks the usual quasi-homologous symmetry. This symmetry breaking is due to the fact 
that A is a fixed parameter. This is an important result, exemplifying the importance of not 
neglecting the cosmological constant. In other words, even being A a quite small constant, its 
presence leads to new information about the vacuum brane system. We conclude emphasizing 
that the quasi-homologous symmetry could, in principle, be restored via an extension of the H 
generator, or by allowing a more elaborate relation between the dark pressure and dark energy. 
These possibilities are currently under investigation. 
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